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We describe neutron scattering experiments performed to investigate the magnetic order and 
dynamics of half-doped Lai.sSro.sCoO-i. This layered perovskite exhibits a near-ideal checkerboard 
pattern of Co 2+ /Co 3+ charge order at temperatures below ~ 800 K. Magnetic correlations are 
observed at temperatures below ~ 60 K but static magnetic order only becomes established at 31 K, 
a temperature at which a kink is observed in the susceptibility. On warming above 31 K we observed 
a change in the magnetic correlations which we attribute either to a spin canting or to a change in the 
proportion of inequivalent magnetic domains. The magnetic excitation spectrum is dominated by 
an intense band extending above a gap of approximately 3 meV up to a maximum energy of 16 meV. 
A weaker band exists in the energy range 20-30 meV. We show that the excitation spectrum is in 
excellent quantitative agreement with the predictions of a spin-wave theory generalized to include 
the full magnetic degrees of freedom of high-spin Co 2+ ions in an axially distorted crystal field, 
coupled by Heisenberg exchange interactions. The magnetic order is found to be stabilized by 
dominant antifcrromagnetic Co 2+ -Co 2+ interactions acting in a straight line through Co 3+ . No 
evidence is found for magnetic scattering from the Co 3+ ions, supporting the view that Co 3+ is in 
the 5 = state in this material. 

PACS numbers: 71.45.Lr, 75.40.Gb, 75.30.Et, 75.30.Fv 



I. INTRODUCTION 

Hole-doped transition-metal oxide antiferromagnets 
exhibit a range of intriguing phenomena, including un- 
conventional superconductivity, metal-insulator transi- 
tions, magnetoelectric behavior and extreme sensitiv- 
ity to external stimuli. Their more unusual properties 
are often found in association with electronically-ordered 
states with nano-scalc periodicity arising from competi- 
tion between electronic, magnetic and structural degrees 
of freedom. As well as being of fundamental interest, 
these highly correlated phases offer the possibility to tune 
the macroscopic properties of a material via selective con- 
trol of its microscopic states. 

At low doping levels, holes introduced into layered an- 
tiferromagnet insulators have a tendency to segregate 
into stripes, as observed for example in layered cuprates 1 , 
nickelates^^ and, most recently, in layered cobaltates 5 . 
Stripes are a form of complex electronic order in which 
the antiferromagnetic structure is modulated by periodic 
arrays of hole-rich antiphase boundaries. With increasing 
doping, stripe-ordered systems can evolve into metallic 
or charge-ordered insulating states depending on the rel- 
ative importance of the electron kinetic energy, Coulomb 
interactions and associated lattice strain. 

At half doping, many transition-metal oxides exhibit 
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an insulating charge-ordered phase, the stability of which 
is generally assisted by cooperative Jahn-Teller distor- 
tions. For layered systems the charge order naturally 
takes the form of a checkerboard pattern. This ordering 
pattern has been reported in the isostructural "214" com- 
ounds La!/2Sr3/ 2 Mn04 (Refjg), La3/ 2 Sr 1 / 2 Co04 (Ref. 
and La3 / / 2 Sr 1 / 2 Ni04 (Ref. Q). Previous investigations 
have revealed a number of interesting features in the or- 
der these compounds. For instance, electronic structure 
calculations for La 1( / 2 Sr3/ 2 Mn04 suggest that there is 
almost no charge separation between the two inequiv- 
alent Mn sites£ In the half-doped nickelate, the checker- 
board charge order rearranges itself spontaneously into a 
stripe-like phase at low temperatures due to spin-charge 
couplingi^iiS Finally, in La3/ 2 Sr!/ 2 Co04, there exists the 
possibility of spin-state transitions^ due to the near de- 
generacy of different terms of the d electron configuration 
in octahedrally-coordinated Co 3+ . In short, these canon- 
ical half-doped perovskites with checkerboard order may 
not be as simple as they first appear. 

In this paper we report the results of neutron scat- 
tering experiments carried out to study the static and 
dynamic magnetic properties of the half-doped cobal- 
tate La3/ 2 Sr!/ 2 Co04. The measurements extend over 
the entire spectrum of cooperative magnetic excitations 
throughout the Brillouin zone, which enables us to char- 
acterize the magnetic excitations fully and to quantify the 
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exchange interactions that stabilize the magnetic ground 
state. The work also sought to establish from the exci- 
tation spectrum whether both of the Co sites are mag- 
netically active, or just one. To this end we carried out 
a rather detailed analysis of the data in terms of a gen- 
eralized spin-wave model that includes the complete set 
of spin and orbital degrees of freedom. We find that the 
measured magnetic excitation spectrum is in excellent 
agreement with the spectrum calculated for the Co 2+ 
site alone. There is no discernible signal in the spectrum 
from the Co 3+ site. 

The crystal structure of Laa^Sr^CoC^ is described 
by the space group IA/mmm, with tetragonal unit cell 
constants a = 3.84 A and c = 12.5 A [see Fig. QTa)]. The 
parent phase La2Co04 is a Mott insulator which orders 
antiferromagnetically below Tn = 275 K. 12 Doping with 
Sr introduces holes into the C0O2 layers, and at half dop- 
ing the 50:50 mixture of Co 2+ and Co 3+ ions crystal- 
lizes at T co ~ 825 K into a checkerboard charge-ordering 
patternfiii 3 - illustrated in Fig. [lib). The signature of 
checkerboard charge order in neutron diffraction data is 
a set of peaks at wavevectors Q co = (h + 0.5, k + 0.5, 1), 
with h, k and I integers (although in the I direction the 
peaks are very broad). Since neutrons do not couple 
directly to charge order, these peaks actually originate 
from associated modulations in the shape of the oxygen 
octahedraiii 3 - 

Magnetic order occurs in Laa^Sr^CoO^ at a temper- 
ature well below the charge ordering temperature. Mag- 
netic susceptibility measurements (Ref [TJ see also Fig. 
U) revealed a broad maximum in the in-plane suscep- 
tibility at ~ 60 K indicative of a build-up of magnetic 
correlations. The data also show a large difference be- 
tween the in-plane (Xab) and out-of-plane (Xe) suscep- 
tibilities, with Xab greater than Xc by at least a factor 
two, revealing strong planar anisotropy. Assuming the 
Co 2+ ions are in the high-spin (HS) state with effective 
spin S = 3/2, as founcU 2 - in La2Co04, Moritomo et alr^ 
concluded from the measured effective moment per Co 
that the Co 3+ ions in La3/2Sr 1 / 2 Co04 are also in the HS 
state (S = 2) and must therefore carry a moment. How- 
ever, this analysis did not take into account the consider- 
able unquenched orbital moment of Co 2+ in the distorted 
octahedral field. Two recent studies, one in which the 
magnetic susceptibility was analysed with a full atomic 
multiplet calculation-^ and the other employing soft x- 
ray absorption spectroscopy to probe the atomic levels 
directly^, have confirmed that the Co 2+ ions are in the 
HS state (S = 3/2) but concluded that the Co 3+ 10ns are 
in the low-spin (LS) state with S = 0, as found for exam- 
ple in LaCoOa. 1 ^ If we accept the weight of evidence in 
favor of the LS state then we can assume that the Co 3+ 
ions are non-magnetic apart perhaps from a small Van 
Vleck moment induced by the exchange field from Co 2+ 
in the ordered phase. 

Neutron diffraction measurements performed by Zal- 
iznyak and coworkers^ confirmed the presence of mag- 
netic order with a gradual onset starting around 60 K. 
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FIG. 1: (Color online) (a) Tetragonal unit cell of 
La2-iSr^Co04. (b) Model for the charge and magnetic order 
within the ab plane of La3/ 2 Sr 1/ / 2 Co04. The shaded rectan- 
gle is the magnetic unit cell. <j> is the angle of the moments 
to the a axis, with positive <j> an anticlockwise rotation. The 
large (yellow) arrow represents the projection of the stacking 
vector t = (1.5, 0.5, 0.5) for the spin-charge order. 



The magnetic order is characterized by a fourfold group 
of magnetic diffraction peaks with slightly incommensu- 
rate wavevectors Q m = (h + 0.5, k + 0.5, ii) ± (0.25 — 
e, 0.25 - e, 0) and (h + 0.5, k + 0.5, h) ± (-0.25 + e, 0.25 - 
e, 0), where h and k are integers, l\ (l 2 ) is an odd (even) 
integer, and e = 0.005 — 0.008 depending on sample 
preparation*^ If the small incommensurability e is ne- 
glected then the simplest magnetic structure consistent 
with the diffraction data would be a collincar antifcr- 
romagnet with ordered moments on the Co 2+ ions and 
propagation vector (0.25, 0.25) within the ab plane, as 
shown in Fig. [ljb). The in-plane orientation of the 
moments (</>) will be addressed later in this paper. It 
has been proposed that the observed incommensurabil- 
ity is caused by stacking faultsJ^ The ideal ordering 
pattern in Fig. [ljb) gives rise to magnetic Bragg peaks 
at in-plane wavevectors (h + 0.5, k + 0.5) ± (0.25,0.25). 
The observed fourfold pattern of peaks arises because 
in tetragonal symmetry a 90deg rotation generates an 
equivalent magnetic structure with propagation vector 
(—0.25,0.25), and in a real sample both wavevector do- 
mains are expected to be present in equal proportion. 

A preliminary report of a subset of the data presented 
here was given in Ref. 1 19L and as far as we are aware 
no other measurements of the magnetic excitations in 
La3/2Sr!/2Co04 have been published. Recently, however, 
there appeared a brief account of the magnetic excita- 
tion spectrum of La3/ 2 Ca 1 / 2 Co04i^ The data indicate 
that the energy scale of the magnetic excitations in the 
Sr-doped and Ca-doped cobaltates is very similar, but 
an important difference is that in Laa^Cai^CoC^ there 
are ordered magnetic moments on both Co 2+ the Co 3+ 
sites M. This increases the number of magnetic modes 
in the excitation spectrum and makes it harder to an- 
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alyze. We also mention that the magnetic excitations 
have been studied in the isostructural half-doped com- 
pounds La 3 / 2 Sr 1 / 2 Ni04 (Ref. [H) and Lai/ 2 Sr 3 / 2 Mn04 
(Ref. H). 

The paper is organized as follows. After giving de- 
tails of the experimental methodology and sample char- 
acterization we describe the results of polarized neutron 
diffraction measurements designed to refine the magnetic 
structure of La 3 / 2 Sr!/ 2 Co04. We then present our inelas- 
tic neutron scattering measurements which represent the 
main body of the work. This is followed by an analysis 
of the magnetic spectrum in terms of a generalized spin- 
wave model. The paper ends with a discussion of the 
results and a summary of the main conclusions from the 
work. 



II. EXPERIMENTAL DETAILS 

Single crystals of La 3 / 2 Sr!/ 2 Co04 were grown in Ox- 
ford by the optical floating-zone method in a four-mirror 
image furnace. The growth took place in a mixed Ar/0 2 
atmosphere at a pressure of 7-9 bar. The feed and seed 
rods were scanned at a rate of 3-4mmhr _1 and counter- 
rotated at 30 r. p.m. Details of the preparation and crys- 
tal growth procedure are very similar to those used for 
Lai_2;Sr x Co0 3 +5 described in Ref. Two large crys- 
tals were prepared for the experiments. These were in the 
form of rods, approximately 10 mm in diameter and up 
to 80 mm in length. Several smaller crystals prepared un- 
der the same conditions were ground up into a fine pow- 
der and used for structural analysis by neutron powder 
diffraction. The general materials (GEM) diffractome- 
ter at the ISIS Facility was used to obtain the powder 
diffraction data. Basic information on the magnetic re- 
sponse of the crystals was obtained from magnetization 
measurements, which were performed on a superconduct- 
ing quantum interference device (SQUID) magnetometer 
(Quantum Design). 

Unpolarized neutron scattering measurements were 
made on the MAPS spectrometer at the ISIS Facility. 
MAPS is a time-of-flight chopper spectrometer with a 
large bank of pixelated detectors. The energy of the in- 
cident neutrons is selected by a Fermi chopper, which 
transmits pulses of neutrons with an energy bandwidth 
of typically 3-5%. All the data presented here were ob- 
tained with an incident energy E\ of 50 meV and a chop- 
per frequency of 350 Hz, giving an elastic energy resolu- 
tion of 2meV. The crystal used on MAPS was a single 
rod with a total mass of 35.5 g. Of this mass an estimated 
24 g was in the neutron beam (which was smaller than the 
size of the crystal). The crystal was mounted with the 
rod axis approximately vertical, and aligned such that the 
tetragonal c axis was parallel to the incident beam direc- 
tion. Spectra were recorded at temperatures of 10 K, 60 K 
and 300 K. Measurements of a standard vanadium sam- 
ple were used to normalize the spectra and place them 
on an absolute intensity scale. 
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FIG. 2: (Color online) Part of the Rietveld refinement for 
tetragonal Lai.sSro.5Co04+i at 60 K. Circles are data mea- 
sured in the backscattering detector bank {{26) = 154.5°) on 
GEM at ISIS. The solid (red) line shows the calculated pro- 
file fit. Tick marks show the positions of allowed reflections, 
and the solid (blue) line below the data shows the difference 
between the data and fit. The background fitted in the re- 
finement was subtracted prior to plotting. The fit parameters 
are given in Table U 



Polarized neutron scattering measurements were made 
on the IN20 and IN22 triple-axis spectrometers at the 
Institut Laue-Langevin. On IN20, two crystals of masses 
6.5 g and 5.5 g cut from the same rod were co-aligned 
with their rod axes parallel. Two settings of the sample 
were used, giving access to the {hkO) and {hhl) planes 
in reciprocal space [we refer crystallographic notation to 
the tetragonal unit cell shown in Fig.flja)]. For the IN22 
experiment we used the 6.5 g crystal on its own and mea- 
sured in just the second orientation. On both IN20 and 
IN22 we used the (111) Bragg reflection of Heusler alloy 
as both monochromator and analyzer and worked with a 
fixed final energy Ef of 14.7 meV. A graphite filter was 
placed after the sample to suppress higher harmonics in 
the scattered beam. Measurements employed uniaxial 
polarization analysis,™ and were made with three config- 
urations of the neutron spin polarization P relative to the 
scattering vector Q: (a) P || Q, (b) P 1 Q with P in the 
scattering plane, and (c) P 1 Q with P perpendicular 
to the scattering plane. Data were collected in both the 
spin-flip (SF) and non-spin-flip (NSF) channels. Details 
of how the six measurements {(a), (b), (c)} x {SF,NSF} 
can be used to separate the components of a magnetic 
structure are given in Ref. [Tol . 



III. RESULTS 

A. Structural characterization 

Neutron powder diffraction data collected on GEM at 
several temperatures were analyzed by Rietveld refine- 
ment using the GSAS suite of programs. 26 Data from 
the three detector banks at highest scattering-angle were 
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refined simultaneously. Refinements were made in the 
space group IA/mmm for the parent structure whose unit 
cell is shown in Fig.[TJa). Attempts to refine the structure 
in the Fmmm space group to allow the possibility of an 
orthorhombic distortion did not lead to any statistically- 
significant improvements in the quality of the fits. In ad- 
dition to the atomic positions shown in Fig.[TJa), we also 
allowed for the possibility of interstitial oxygen at site 
(0.5,0,0.25), where excess oxygen has previously been 
founds in La 2 Co04 + a. In order to achieve convergence 
it was necessary to fix the thermal parameter (J7 iso ) of 
this third oxygen site to be the same as that of the sec- 
ond oxygen site. 

As an example, Fig. [5] shows the neutron powder 
diffraction pattern and refined profile fit for the data 
recorded in the highest-angle detector bank with the sam- 
ple at a temperature of 60 K. The parameters for the re- 
finement are listed in Table U together with those for the 
other temperatures at which measurements were made. 
The data at 2K will have contained magnetic Bragg 
peaks which were not allowed for in the refinement. This 
probably explains the slightly higher value of R wp for this 
temperature. Values for the lattice constants and atomic 
positions are in good agreement with previous data^ 
To within experimental accuracy the refinements show 
no deviation from stoichiometry in the oxygen content 
(i.e. (5 = 0), and no evidence for any deviation from the 
nominal La:Sr ratio of 3:1. 



B. Magnetic structure 

We first review the main characteristics of the mag- 
netic and charge order in Lai.sSro.sCoC^ with the aid 
of the polarized-neutron diffraction measurements pre- 
sented in Fig. [3] All the scans shown in this figure con- 
tain raw data recorded at 2 K on IN20 with the neutron 
polarization parallel to the scattering vector (P || Q). 
In this configuration the spin-flip scattering is entirely 
magnetic and the non-spin-flip scattering is entirely non- 
magnetic (i.e. structural). Since the neutron polarization 
is not perfect (the flipping ratio was ~ 20) there is a small 
amount of leakage from one channel to the other. 

Figure 0(a) is a diagram of part of the (h, k, 0) plane 
in reciprocal space showing the in-plane wavevectors cor- 
responding to the magnetic and charge order. Repre- 
sentative scans along various in-plane and out-of-plane 
directions are shown in Figs. OHb)-(e). These data 
confirm that the charge order scattering is peaked at 
wavevectors Q co = (h + 0.5, k + 0.5, 1), with h, k 
and I integers, and that the magnetic order is char- 
acterized by four slightly-incommensurate wavevectors 
Q m = (h + 0.5, k + 0.5, 1) ± (0.25 - e, 0.25 - e, 0) and 
(h + 0.5, k + 0.5,0 ± (-0.25 + e,0.25 - e,0), where h 
and k are integers and I is an odd integer for the first 
pair of Q m and an even integer for the second pair. In 
our crystal e = 0.005 ± 0.001. The peaks are consider- 
ably sharper in scans parallel to the ab plane [e.g. Figs. 



|3jc) and (d)] than in scans made in the out-of-plane 
direction [Figs. [3fb) and (e)], but in all directions the 
peak widths are broader than the resolution. The cor- 
relation lengths for the charge and magnetic orders esti- 
mated from £ = 1/r, where T is the half- width at half- 
maximum, are = 23±2 A, £r Q = 8±1 A, ^ b = 52±2 A 
and £^ = 12 ±1 A. These values are consistent with those 
reported by Zaliznyak et al. 7 apart from which we 
find to be ~ 30% smaller. This could be because we did 
not attempt to correct the peak width for experimental 
resolution. 

The two in-plane scans shown in Fig. [3Jc) are at I = 0, 
which is a minimum of the intensity modulation along 
(0, 0, 1) — see Fig. [21(e). The widths of the peaks in 
Fig. [2fc) are roughly twice those in Fig. EKd) which are 
at a maximum of the (0, 0, 1) intensity modulation. This 
suggests that in addition to the dominant magnetic order 
there also exists a small volume fraction of magnetically 
ordered regions with an average size of ~ 25 A in the ab 
plane which are not correlated to the magnetic order on 
the adjacent layers. 

It is reasonable to assume that the short correlation 
lengths in the c direction are caused by the existence 
of different stacking sequences with similar energy. Since 
both the charge and magnetic diffraction peaks are found 
at integer I the majority stacking is periodic in the lattice. 
However, there are several ways in which the structure 
on adjacent layers can be related. To distinguish these 
experimentally let us consider the case where the mag- 
netic structure is collinear, so that the spin-charge order 
on one layer (say z = 0) can be related to that on the 
adjacent layer (z = ^) by a translation t. The structure 
factor for the spin and charge diffraction peaks then con- 
tains a factor l+exp(iQ-t) = l+exp{2wi(ht x +kt y +lt z )}, 
where t x ,t y ,t z are the components of t written as frac- 
tional coordinates along the crystallographic a, b and c 
axes. 

We consider first the case t = (—0.5,0.5,0.5). For 
the structure shown in Fig. |TJt>) this stacking leads to 
systematic absences at magnetic ordering wavevectors of 
the type (0.25,0.25,0, (0.75,0.75,0, etc, when I is odd 
and for magnetic wavevectors of the type (0.25,0.75,0, 
(0.75,0.25,/), etc [associated with the twin obtained by 
rotating the structure in Fig. QJb) by 90 deg — see Fig. 
131(a)], when I is even. These predictions are inconsis- 
tent with experiment — see Fig. |3Je) . We next consider 
t = (0.5,0.5,0.5). In this case there are no magnetic ab- 
sences, again inconsistent with experiment. We finally 
consider t = (1.5,0.5,0.5). For this case the systematic 
absences are at (0.25,0.25,0, etc, when I is even, and 
at (0.25,0.75,0, etc, when I is odd. This is consistent 
with the observations. The twinning ensures that there 
are no absences in the charge-order peaks, again consis- 
tent with experiment^ We conclude, therefore, that the 
most likely stacking vector is t = (1.5,0.5,0.5) (or its 
equivalent). This vector is shown in Fig.QJb). With this 
stacking, spins in the z = | layer are antiparallel to the 
closest spins in the z = and z = 1 layers. 
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TABLE I: Structural parameters for Lai.sSro.sCoO^i at temperatures between 2K and 300 K refined from neutron powder 
diffraction data. The refinements were performed in the tetragonal space group I4/mmm of the parent structure with atomic 
positions La/Sr 4e (0,0, z La ), Co 2a (0,0,0), O(f) 4c (0.5,0,0), 0(2) 4e (0,0, 202), 0(3) 4d (0.5,0,0.25) in lattice units. The 
0(3) position is an interstitial site, n is the occupancy of each site per formula unit and Ui so is the isotropic temperature factor. 
The occupancy of the Co site was not refined (n = f), and La/Sr site was constrained such that iiLa + ns r = 2. The numbers 
in parentheses are statistical errors in the last digit of the refined parameters. R wp is the weighted profile residual function. 
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The fact that the magnetic and charge correlation 
lengths along the c axis are so short suggests that the 
majority stacking is only slightly more favorable ener- 
getically than other possible stacking sequences. There- 
fore, the propagation of the structure along the c axis is 
presumably interrupted frequently by stacking faults in 
which adjacent layers are related by different t vectors. 

Let us now turn to the temperature evolution of the 
magnetic order. In Fig. 3] we show the magnetization of 
Lai.5Sro.5Co04 as a function of temperature, recorded 
with a relatively low measuring field H of 100 Oe. The 
data are in very good quantitative agreement with the 
measurements reported in Refs. ITT1 and [Til both of which 
employed a significantly higher measuring field. The 
magnetization exhibits strong XK-like anisotropy, as 
shown in the inset of Fig. 31 and there is a broad max- 
imum in the in-plane response centered at about 60 K. 
At lower temperatures there is a splitting between the 
field-cooled (FC) and zero-field-cooled (ZFC) data. The 
splitting is largest for H \\ ab. 

Our measurements, however, reveal an additional fea- 
ture. At approximately 31 K there is a sharp kink in 
the data. This kink marks the temperature below which 
the FC-ZFC splitting begins to open up most rapidly. 
The observation of this kink prompted us to perform 
polarized-neutron diffraction measurements to investi- 



gate whether there might be a change in the magnetic 
structure at 31 K. 

We followed the approach described in Ref. [llj (see 
also Ref. |28[ ). SF and NSF intensities were recorded at 
two magnetic Bragg peaks using three orthogonal direc- 
tions of the neutron polarization P. Measurements were 
made on both IN20 and IN22, and corrections were ap- 
plied to the measured intensities to compensate for the 
imperfect polarization. The two Bragg peaks used were 
Qi = (0.25 + e, 0.25 + e, 7) and Q 2 = (1.25+e, 1.25 + e, 1). 
These were chosen because they make small angles (less 
than lOdeg) to the c axis and [110] direction, respec- 
tively, which reduces the uncertainty in the values of the 
spin components derived by this method. For example, if 
Qi could be chosen exactly parallel to c then the Bragg 
peak intensities measured at Qi would be independent 
of the magnetic component along the c axis. 

Because Qi and Q 2 lie in the (h 7 h, I) plane the mag- 
netic scattering at these wavevectors is naturally de- 
scribed in terms of the intensities scattered by the pro- 
jection of the ordered moments along the orthogonal di- 
rections [1,1,0], [1,-1,0] and [0,0,1]. We call these in- 
tensities /no, inn and J c , respectively. The expressions 
111 Table I of Ref. 10 can be used to determine the ratios 
^110/^110 an d Ic/{hw + I110) from the sets of measure- 
ments at Qi and Q 2 . The latter ratio was found to be less 
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FIG. 3: (Color online) Polarized neutron diffraction from 
magnetic and charge orders in Lai.sSro.sCoOzt. (a) Dia- 
gram of the (h, k, 0) reciprocal lattice plane showing the in- 
plane wavevectors of the magnetic order (two twins, unfilled 
and filled circles marked B and C, respectively) and charge 
order (crossed circle marked A). The arrows indicate scan 
directions, (b) /-dependence of the charge-order peak at 
(0.5,0.5,7). (c) and (d) In-plane scans through the mag- 
netic ordering wavevectors in the directions indicated, (e) 
/-dependence of the intensity at the magnetic wavevector 
marked B. Shaded grey regions contain points contaminated 
with scattering from aluminium. All scans shown in this fig- 
ure were made on IN20 with P | Q and recorded in both the 
spin-flip (SF, filled blue circles) and non-spin-flip (NSF, open 
red squares) channels. No corrections have been applied for 
the imperfect neutron polarization. The sample temperature 
was 2K. 
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FIG. 4: (Color online) Temperature dependence of the mag- 
netization of Lai.5Sr .5CoO4. The main frame shows mea- 
surements made with the applied field H of strength 100 Oe 
(fioH = 0.01 T) applied parallel to the ab plane. Red filled cir- 
cles and blue open squares show data from field-cooled (FC) 
and zero-field-cooled (ZFC) measurements, respectively. In- 
set: data up to T — 350 K showing measurements with both 
H II ab and H || c. 
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FIG. 5: (Color online) Development of magnetic order in 
Lai.5Sro.5Co04. (a) Ratio of the intensities scattered by the 
projection of the ordered moments along the [110] and [110] 
directions. Data are from IN20 (filled red circles) and IN22 
(filled blue squares), (b) Intensity of the (0.75 — e, 0.75 — e, 1) 
magnetic Bragg peak. The broken grey line marks the posi- 
tion of the kink in the magnetization (see Fig. [4}. 



than 0.01 at all temperatures. The assumption that the 
intensities are proportional to the squares of the ordered 
moments constrains the angle of the moments to the ab 
plane to < 5 deg. Given the strong planar anisotropy it 
is safe to assume that the moments lie in the ab plane. 

In Fig.[5£a) we plot the ratio Iuq/Iiiq as a function of 
temperature. At temperatures below 30 K, Tno/iiio is 
approximately constant with a value of 0.4. On warming 



above 30 K this ratio gradually increases until at ~ 50 K 
it approaches 1. This indicates one of two things. Either 
(i) there are two (or more) spin domains which contribute 
to each magnetic Bragg peak, and on cooling one of these 
domains becomes preferentially populated, or (ii) there 
is a spin-canting transition reminiscent of that proposed 
to occur in the isostructural layered nickelates i 10 i 29 ' 30 i 31 
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FIG. 6: (Color online) Spin excitation spectrum of 
Lai.sSro.sCoCU. (a) Map of the intensity within the (h,k) 
plane averaged over the energy range 7 — llmeV, measured 
with a sample temperature of T = 10 K . The circles cen- 
tered on the magnetic ordering wavevectors are due to scat- 
tering from dispersive magnetic excitations, (b) Energy-Q 
slice showing the dispersion along the line Y in (a), (c) Energy 
scan at the magnetic ordering wavevector (0.75 — e, 0.25 + e, 0) 
with neutron polarization analysis to separate the magnetic 
scattering in the spin-flip (SF) channel from the non-magnetic 
scattering in the non-spin-flip (NSF) channel. The data were 
obtained on IN20 at T = 2K and reveal an energy gap of 
~ 3meV in the magnetic spectrum, (d) Scans along the 
line X in (a) at temperatures of 10 K and 300 K averaged 
over the energy range 28 — 32 meV confirming the existence 
of magnetic modes at this energy. The data in (a), (b) and 
(d) were recorded on MAPS, and the intensity is in units of 
mb sr _1 meV -1 f.u. -1 , where "f.u." stands for "formula unit" 
(of Lai.5Sr . 5 CoO4). 



We consider these possibilities further in Section fVl 



Figure [5^b) shows the temperature dependence of the 
intensity of the (0.75 — e, 0.75 — e, 1) magnetic peak. The 
slow increase in intensity on cooling starts in the vicinity 
of the broad hump in the magnetization at ~ 60 K, con- 
firming that the hump is associated with the build-up of 
magnetic correlations. We have indicated on Figs. 02a) 
and (b) the temperature at which the kink is observed 
in the magnetization. There is no obvious anomaly in 
the magnetic peak intensity at this temperature, but the 
data in Fig. Ela) indicate that the 31 K kink marks the 
temperature below which the magnetic structure stops 
changing. 



C. Magnetic excitations 

Figure [6] provides an overview of the magnetic spec- 
trum. The intensity map presented in Fig. E3a) is a slice 
from a MAPS data volume, averaged over the energy 
range 7— llmeV and projected onto the (h, k) plane. 
The image shows dispersive magnetic excitations emerg- 
ing from the magnetic ordering wavevectors. The rings of 
scattering correspond to the intersection of the constant- 
energy slice plane with the spin-wave cones. Figure 
|H{b) is an energy- wavevector slice through the same data 
volume showing the dispersion along a line parallel to 
(h, h, 0). The spectrum is dominated by an intense band 
of scattering extending up to 16meV. There is also a 
weaker band of scattering in the range 20 — 30meV which 
is more diffuse than that of the lower band but which dis- 
perses with the same period. 

Figure [He) is an energy scan performed on IN20 with 
the wavevector fixed at the magnetic ordering wavevec- 
tor Q m — (0.75 — e, 0.25 + e, 0). The two sets of points 
are the neutron spin-flip (SF) and non-spin-flip (NSF) 
signals. In the P || Q configuration used for this mea- 
surement the SF channel contains only magnetic scat- 
tering, so the large signal in the SF data is scattering 
from magnetic excitations and reveals a gap of approx- 
imately 3 meV in the low-energy spin-wave band. Hav- 
ing observed the gap we performed constant-energy scans 
(not shown) along (0.75 — e,0.75 — e,l) from I — to 
2 at energies of 2meV and 4meV, i.e. just below and 
just above the gap. At 2 meV there remained a mod- 
ulation in the magnetic scattering along the scan with 
a broad maximum at the magnetic Bragg peak position 
(0.75 — e, 0.75 — e, 1), whereas at 4meV the scan was fea- 
tureless to within the experimental precision of ~ 5%. 
This implies that for energies above 4 meV the magnetic 
dynamics are completely uncorrelated along the c axis 
and the spectrum can be considered as two-dimensional. 

To help understand the origin of the gap we per- 
formed neutron polarization analysis at energies of 2 meV 
and 4meV. For each energy we recorded the signal in 
the same three polarization channels and at the same 
two wavevectors Qi = (0.25 + e,0.25 + e, 7) and Q 2 = 
(1.25 + e, 1.25 + e, 1) as used to analyze the magnetic or- 
der (see Section Till Bj) . The signal in each polarization 
channel can be written in terms of the response func- 
tions S^CQj.cd), S^iQ^u) and S zz (Qj,uj) for mag- 
netic fluctuation components parallel to [110], [110], and 
parallel to the c axis [z direction), respectively. Here, u> 
is the neutron energy transfer and j = 1,2 indexes the 
two magnetic wavevectors. Applying the same analysis 
as used to separate diffraction from different components 
of the ordered moments 10-28 we find at T = 2K, 



S zz (Qj,2 meV) /S 110 (Q rf , 2 meV) 
S zz (Qj , 4 meV) /S ll ° (Qj , 4 meV) 
S 1W (Qj , 2 meV) /S ll ° (Qj,2 meV) 
S 1W (Qj , 4 meV) /S* 110 (Qj , 4 meV) 



-0.04 ±0.06 
-0.02 ±0.04 
1.19 ±0.07 
1.00 ±0.05 . 
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FIG. 7: (Color online) Intensity maps of measured and simulated neutron scattering from Lai.sSro.sCoGi. The data are from 
MAPS, and the simulations are obtained from the many-level spin-wave model discussed in the text with Co 2+ — Co 2+ exchange 
interactions J — 1.4 meV and Ji = J2 = [see Fig. H^b)]. ( a )~(f) Constant-energy slices averaged over the energy ranges 
indicated, (g)-(l) Energy-Q slices showing the dispersion along the three symmetry directions indicated. The intensity scale 
is in units of mbsr -1 meV -1 f.u. -1 . 



The first two ratios show that the spin fluctuations at 
both 2meV and 4meV are restricted to the ab plane. 
This is consistent with the strong XF-like anisotropy of 
this system and implies that the gap is due to a small 
single-ion or exchange anisotropy within the ab plane. 
The third and fourth ratios show that the strengths of the 
fluctuations in the [110] and [110] directions are roughly 
equal in this energy range. 

In Fig. [S^d) we show cuts through the MAPS data av- 
eraged over the energy range 28meV to 32meV, which 
is near the top of the weak upper band of excitations. 
The scan at 10 K shows a sinusoidal intensity modulation 
along the cut direction [line X in Fig. [5(a)], whereas the 
scan at 300 K is flat and always below the 10 K data. The 
observation that the signal goes away on raising the tem- 
perature from 10 K to 300 K shows conclusively that it is 
magnetic in origin, because the scattering from phonons 
would increase in intensity. 

To give a more comprehensive picture of the magnetic 
excitation spectrum we present in Fig. [7] a further se- 
ries of slices through the MAPS data. The data have 
been averaged over symmetry-equivalent directions to 



improve the statistics, and are plotted with respect to 
the two-dimensional reciprocal lattice of the C0O2 layers 
indexed by (h,k). In the time-of-flight method the out- 
of-plane wavevector component / varies with energy and 
also varies across the detector for a fixed energy. How- 
ever, because the magnetic dynamics is two-dimensional 
(for energies greater than ^4 meV) the magnetic disper- 
sion does not depend on I and the intensity of the spec- 
trum varies only slowly with I. The I dependence of the 
intensities of the modes is correctly taken into account in 
the models presented in Section ITVl 

Figures 0a)-(c) are constant-energy slices averaged 
over a 2meV interval centered on 5meV, 9meV and 
13meV. These illustrate how the lower spin- wave band 
disperses away from the magnetic zone centers. Fig- 
ures|7(g)-(i) display three energy-wave vector slices taken 
along different symmetry directions in the Brillouin zone. 
The intense lower band and weak upper band are clearly 
seen in each of these slices. 

To facilitate comparison with models we took cuts 
through the MAPS data set along several symmetry di- 
rections and extracted points describing the dispersion. 
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The magnetic dispersion obtained from the MAPS data 
is displayed in the upper panel of Fig. \E\ The points are 
obtained either from wavevector cuts at fixed energy, or 
from energy cuts at fixed wavevector. Peaks in the cuts 
were fitted with Lorentzian functions on a linear back- 
ground. The points representing the energy gap at the 
magnetic zone centres were estimated from the data in 
Fig. He). 



IV. ANALYSIS 

The essential physics governing the magnetic proper- 
ties of Lai.5Sr .5CoO4 can be understood with reference 
to Fig. [9j which is an energy-level diagram for the high- 
spin state of Co 2+ (3d 7 , L = 3, S = 3/2) in an axially- 
distorted octahedral crystal field. A perfect octahedral 
field splits the orbital levels into two triplets and a sin- 
glet, and the small axial elongation of the octahedron 
splits each triplet into a singlet and a higher-lying dou- 
blet. Inclusion of spin-orbit coupling lifts the fourfold 
spin degeneracy of each of the orbital levels, splitting 
the ground state into two doublets. For realistic crys- 
tal field and spin-orbit interactions the two doublets are 
separated by approximately 25 meV, as we show later. 
Because the orbital ground state in the main octahedral 
field is a quasi-triplet it contains significant unquenched 
orbital angular momentum which is responsible for the 
strong planar anisotropy. As we shall see, the magnetic 
spectrum in the energy range probed in this study in- 
volves the lowest two doublets and is strongly influenced 
by the in-plane anisotropy and the unquenched orbital 
angular momentum. 

The neutron scattering cross section for spin-only scat- 
tering, or for spin and orbital scattering in the dipole 
approximation, may be written^ 



a/3 



Q a Qp)S aP (Q t u>) t (1) 



where (jro/2) 2 = 72.8 mb, f(Q) is the magnetic form fac- 
tor of Co 2+ , exp(— 2W) is the Debye- Waller factor which 
is close to unity at low temperatures, Q a is the a compo- 
nent of a unit vector in the direction of Q, and 5 a ^(Q, u) 
is the response function describing a/3 magnetic correla- 
tions. In calculating the cross-section we average over an 
assumed 50:50 ratio of the two magnetic domains related 
to one another by a 90deg rotation. 

Putting the orbital angular momentum to one side for a 
moment, we can as a first approximation follow the stan- 
dard route and attempt to describe the low-energy mag- 
netic dynamics in the antiferromagnetic phase in terms of 
linear spin-wave excitations of an effective spin-^ model 
for the ground state doublet. We assume the collinear 
magnetic structure shown in Fig. [1] with spins along the 
a axis (0 = 0). We consider a spin-only Hamiltonian and 




0.10.20.30 4 30 4 0.5 0.4 0.3 0.2 0.1 D.1 0.2 0.3 4 O S 
(0.25. 1 (F), 0.5) (h. h) (h,0) (r.l.u.) 



p, 
■i 

2 

v 

\ 4 
' 2 
' 
100 

50 



if 








« , 
















Si 








A 
\ 

f \ 

22 




x 





0.1020304 
(0.25, k) 



030405 04 03 02 0.1 D102030405 
(h, 0.5) (h. h) (h.0) <i.l.u.) 



FIG. 8: (Color online) Dispersion of magnetic excitations in 
Lai.sSro.sCoCU. The diagram at the top shows the path in 
reciprocal space along which the dispersion is plotted. The 
middle figure shows the measured and calculated dispersion. 
Symbols are from fits to the experimental data. Circles are 
from constant-energy cuts, squares and diamonds are from 
constant-wavevector cuts. The lines are calculated from the 
many- level spin-orbital model with Co 2+ — Co 2+ exchange in- 
teractions J = 1.4 meV and Ji = J2 = [see Fig. H]T)]. 
The lower figure displays the calculated response functions 
S aa (Q,Lj) for each mode shown in the dispersion plot. The 
labels on the curves give the aa components for the case where 
the ordered moments point along x. 



incorporate the magnetic anisotropy via anisotropic ex- 
change interactions between effective S = h spins. The 
Hamiltonian may be written 



nt \ ^ \ ^ roc na Qa 

n — 2-^i i k 3 k ' 

(Jk) a 



(2) 



The first summation is over Co 2+ — Co 2+ pairs with each 
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FIG. 9: Single-ion energy levels of the F term of Co + in 
Lai.5Sro.5Co04, showing how the term splits when succes- 
sively smaller interactions are added. The splitting caused by 
the axially-distorted crystal field is obtained with coefficients 
B% = 13meV, B° A = -1.4 meV and B% = -8meV for the 
Stevens operator equivalents 0\ , 0% and 0% . These values are 
approximately twice those calculated from the point-charge 
model. The right-most scheme includes spin-orbit coupling 
and shows the splitting of the lowest orbital quasi-triplet on 
an expanded vertical scale. The spin-orbit coupling parame- 
ter is A = — 18.7meV. 
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FIG. 10: (Color online) Magnetic dispersion parallel to the 
(ft., 0) direction in Lai.sSro.sCoO^ The symbols are points 
determined from energy cuts at constant wavevector through 
the MAPS data. Full and broken lines are in-plane and out- 
of-plane modes of the effective spin-1 ground-state doublet 
calculated from Eq. ([2]) by linear spin-wave theory with the 
following exchange parameters: (i) SJ = 3.68 meV, J\ = Ji = 
-J/2; (ii) SJ = 3.23 meV, SJi = SJ 2 = -0.15 meV; (hi) 
SJ = 2.85 meV, Ji = J2 = J/2. The parameters in (ii) have 
been adjusted to give the best fit to the data. In each case the 
exchange-anisotropy parameters were e = 0.03 and S = 0.65. 



pair counted only once, and the second summation is 
over the spin components a = x,y,z. The J" k are the 
exchange parameters. We include only J, J\ and J2 as 
defined in Fig. [JJ J acts in a straight line through the 
Co 3+ site, whereas J\ and J2 have 90deg paths. As we 
shall see, J\ and J2 are very much smaller than J, and so 
as a simplification we take J\ and J2 to be isotropic. For 
later convenience we write the anisotropy in J in the form 
J* = J(l + E ), jv = J and J z = J(l - 6) so that e and 
8 parameterize the degree of in-plane and out-of-plane 
anisotropy, respectively. 

The theoretical expressions for the spin-wave disper- 
sion and response functions are given in the Appendix, 
Eqs. ([I])-©. For each wavevector there are two modes 
which are degenerate in the absence of anisotropy. Inclu- 
sion of easy-plane anisotropy lifts the degeneracy so that 
one mode (oji) is associated with in-plane fluctuations 
and the other mode (0^2) is associated with out-of-plane 
fluctuations. At the magnetic Bragg peak wavevectors 
(0.25,0.25), (0.75,0.75), etc. the out-of-plane mode is 
gapped while the in-plane mode is gapless. The addition 
of a small in-plane anisotropy gaps the in-plane mode too. 
In our case we require a large out-of-plane anisotropy to 
describe the data. To see this we refer to the scattering 
intensity maps, e.g. Fig. [7^g), and the dispersion data 
in Fig. [H] If the anisotropy were small then the disper- 
sion would reach a maximum at, or very close to, the 
magnetic zone boundaries (0.125,0.125), (0.375,0.375), 
etc. However, the measured dispersion reaches a maxi- 
mum not at these positions but at the zone centres (0, 0), 
(0.5,0.5), etc. This behavior can be reproduced in the 
linear spin-wave model only if there is substantial easy- 



plane anisotropy, i.e. if J z is significantly less than J x 
and jv. 

We found that the observed dispersion for the lower 
energy modes can be described very well throughout the 
Brillouin zone by the effective spin-i spin-wave model 
with parameters SJ = 3.3 meV, J\ — J2 = 0, 5 = 0.65 
and e — 0.03. The zero values of the 90deg exchange 
parameters J\ and J2 means that the magnetic struc- 
ture consists of two uncoupled, interpenetrating, square- 
lattice antiferromagnets with spacing 2a on the C0O2 
layers. This is surprising and has important implications 
for the stability of the magnetic structure, so let us ex- 
amine the evidence for this finding carefully. 

First of all we show in Fig. [TU] the measured dispersion 
along the (h, 0) line in reciprocal space. If the tetragonal 
symmetry is preserved in the charge-ordered phase then 
one expects J\ — J2, which will tend to frustrate the 
magnetic order. If J\ and J 2 are both antiferromagnetic 
then the effect on the dispersion relation is to soften the 
in-plane mode at the zone boundaries along the (h, 0) and 
(0, k) directions. If J\ and J2 are both ferromagnetic then 
the softening is instead at the reciprocal lattice vectors. 
In both cases, therefore, we should observe a difference 
between the magnon energy at (0.5,0) and at (1,0). It 
can be seen from the data points in Fig. [TU] that if such a 
difference exists then it is very small. To be more quanti- 
tative we plot on Fig [10] the dispersion curves calculated 
from linear spin-wave theory for the cases J\ = J 2 = J/2 
(J, Ji, J2 all antiferromagnetic) and J\ — J2 = — J/2 (J 
antiferromagnetic, Ji, Ji ferromagnetic), adjusting the 
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value of J to accord with the data at (1,0). These curves 
fall respectively well below and well above the experi- 
mental data, showing that the magnitudes of J\ and Ji 
must be considerably less than half that of J. The third 
curve plotted in Fig [10] is obtained with the parameters 
SJ = 3.23 meV, SJi = SJ 2 = -0.15 mcV, and with 
<5 = 0.65 and e = 0.03 as found earlier. This set of 
parameters gives the best fit to the data given the con- 
straints, but is only a marginal improvement on the case 

J t = J 2 = 0. 

In the more general case where charge-order leads to 
orthorhombic symmetry, J\ and Ji could be different. 
Let us therefore consider the case where J and J\ are an- 
tiferromagnetic and J 2 is ferromagnetic, consistent with 
the structure in Fig. [TJb) . This also creates differences 
between the magnon energies at (0.5, 0) and at (1, 0), and 
in addition it increases the energy of the in-plane mode at 
wavevectors such as (0.25,0.75) and (0.75,0.25). Neither 
of these effects is observed experimentally. 

The spin-wave model just described succeeds in giv- 
ing a good description of the dispersion of the lowest 
energy branches of the magnetic excitation spectrum 
and provides compelling evidence that the magnetism in 
Lai.5Sr .5CoO4 is dominated by 180 deg Co 2+ — Co 2+ ex- 
change interactions. However, this model fails in two 
respects. First, it does not account for the band of mag- 
netic scattering observed in the range 20 — 30 meV, and 
second it cannot predict the correct intensities of the 
modes because it ignores orbital angular momentum. 

The effective spin- 1 model takes for its basis the eigen- 
functions of S x , where x is the spin quantization direc- 
tion. When single-ion anisotropy is strong and the true 
spin is greater than i, as we have here, there can be sig- 
nificant admixing of the basis states. This means that 
excitations to higher single-ion levels can propagate and 
can be observed by neutron scattering. Moreover, when 
the single-ion states contain a non-negligible orbital com- 
ponent this needs to be included for an accurate calcu- 
lation of the neutron cross section, because the neutron 
couples to both the spin and the orbital angular momen- 
tum. 

To achieve a more complete description of the mag- 
netic spectrum we use a more realistic linear spin- 
wave formulation which takes as its basis the product 
states \L X S X ) and determines the single- ion states self- 
consistently through the action of the crystal and ex- 
change fields and the spin-orbit coupling. The method 
follows closely the approach described in Ref. 33 in con- 
nection with the dispersive magnetic excitations of an- 
tiferromagnetic KC0F3. Similar calculations have been 
done for CoO (Refs. M and HI) and for CoF 2 (Ref. [H), 
the latter though without explicit inclusion of the orbital 
angular momentum. 

The Hamiltonian for the coupled 3d states of Co 2+ 
(3d 7 , L — 3, S — 3/2) in the distorted octahedral ligand 



field in Lai.5Sro.5Co04 is taken to be 



H = 22 JjkSj ■ s fe + ^2 



(3) 



The first term is an isotropic Heisenberg superexchange 
interaction between true spins S = §, and the second 
and third terms describe the single-ion crystal (ligand) 
field and spin-orbit interactions, respectively. The crys- 
tal field interaction is assumed to be the source of mag- 
netic anisotropy, and any anisotropy in the exchange 
interactions is neglected. The O™ are Stevens opera- 
tor equivalents with B™ the corresponding crystal-field 
parameters, and A is the spin-orbit coupling parame- 
ter. We used the value A = — 18.7meV recently ob- 
tained from optical spectroscopy of CoOi^I The axially- 
distorted octahedral crystal field is represented by the 
Stevens operators O®, O4 and Of, where the quantiza- 
tion direction is taken along the tetragonal axis. We 
used the point charge models to make a first estimate 
of the crystal-field parameters, using the data in Ta- 
ble U for the average Co-0 bond lengths and the re- 
sults of Ref. [H for the oxygen displacements due to 
charge order, which slightly reduces the axial distortion 
of the octahedron. This predicted a total splitting of the 
3d 7 4 F term of about 800 meV (the small admixture of 
higher terms was neglected). The crystal- field splitting in 
Lai.5Sro.5Co04 has not been measured, but optical mea- 
surements of several compounds containing octahcdrally- 
coordinated CoC>6 complexes show a typical splitting of 
about 2 eV i 39 ' 40 This indicates that the point-charge pa- 
rameters are about a factor 2 too small, so for an ini- 
tial estimate we doubled the point-charge values, giving 



B% = 9.3 meV, Bl 



-1.35 meV and B\ = -8.0 meV. 



We note that with this crystal field the ground state 
has XK-like anisotropy, as observed experimentally. The 
ground to first-excited doublet-doublet splitting is not 
very sensitive to the B® and B\ parameters, but varies 
strongly with B®. We therefore adjusted B^ to obtain a 
good fit to the experimental data. 

The method for diagonalizing ([3]) is described in the 
Appendix. To compare the many-level spin-orbital 
model with the data we fixed J\ = J2 = 0, since we have 
already established that these parameters are very small, 
and allowed only the parameters J and B® to vary. The 
best agreement was obtained with B® = 13.0 meV and 
,/ = 1.4 meV. A small anisotropy- field term H a • S with 
H a directed along the x axis and of magnitude 0.22 me V 
was added to ([3]) to fix the direction of the ordered mo- 
ments and to reproduce the observed in-plane spin gap of 
~3 mcV. The dispersion relations and response functions 
calculated with these parameters are shown in Fig. [5] 
The calculated dispersion for the low-energy branches is 
virtually indistinguishable from that calculated from the 
effective spin-i model described earlier and agrees very 
well with the experimental data. We note that when 
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Ji = J% = the magnetic dispersion is the same for the 
two twins of the assumed magnetic structure [Fig. QJb)] 
whose propagation vectors are related by a 90° rotation. 
This means that the analysis of the excitation spectrum 
is not affected by the twinning. 

The many-level spin-orbital model allows us to go be- 
yond the lowest two magnon branches and examine the 
modes derived from the upper doublet in the single-ion 
spectrum (Fig. [9]). The model predicts two pairs of dis- 
persive bands between 25 meV and 30 meV as shown in 
the middle panel of Fig. [8] In the lower panel of Fig. [8] it 
can be seen that two out of these four modes have non- 
negligible response functions. Interestingly, the mode 
with strongest intensity is a longitudinal excitation, la- 
beled xx in Fig.0 and is therefore very different in char- 
acter to a conventional spin precession wave. 

In Fig. [7] we show neutron scattering intensity maps 
calculated with the many-level spin-orbital model along- 
side the corresponding experimental data. The presented 
intensity is {ki/kj) d 2 cr/dild£'f , i.e. the scattering cross 
section in absolute units as defined in Eq. [T] multiplied 
by the factor ki/kj. The simulations properly take into 
account the variation in the scattering vector Q with po- 
sition on the detector and with energy, and the intensity 
has been averaged over equivalent 90deg magnetic do- 
mains. The simulated spectra were convoluted with a 
Lorentzian broadening function with a full-width at half 
maximum of 4 meV to take into account the spectrometer 
resolution (~ 2 meV) and intrinsic broadening, and have 
been multiplied by a constant scale factor of 0.3 to obtain 
quantitative agreement with the experimental data. We 
estimate that absorption and self-shielding in the sample 
accounts for a factor of about 0.5, and a further reduc- 
tion of order 10% may be expected due to zero-point 
magnetic fluctuations not included in the model. There- 
fore, the known corrections account for a scale factor of 
about 0.45, which is somewhat larger than the applied 
scale factor of 0.3. 

Overall, the simulations from the many-level spin- 
orbital model provides an excellent description of the 
observed magnetic excitation spectrum in the measured 
energy range. With a realistic crystal field and only two 
free parameters in the Hamiltonian (J and i?^) plus an 
overall scale factor the simulations reproduce the disper- 
sion and intensities of the all the observed modes in the 
spectrum extremely well. 



V. DISCUSSION 

The experiments presented here have shown that both 
the magnetic order and the magnetic excitation spec- 
trum of charge-ordered Lai.sSro.sCoCU can be under- 
stood in terms of a square, two-sublattice, collinear an- 
tiferromagnet, with ordered moments localized on the 
Co 2+ ions and pointing in the C0O2 layers. The two 
lowest-frequency branches of the magnetic excitation 
spectrum are conventional spin-precession waves, but 



higher-frequency dispersive modes are of different char- 
acter. In particular, a branch observed at ~ 25 meV 
is found to correspond to propagating magnetic fluctua- 
tions which are longitudinal relative to the ordered mo- 
ment direction. 

To arrive at this point, proper account has had to be 
taken of the orbital component in the ground state, which 
is not negligible. For example, from the self-consistent 
solution of the mean-field equations ||5J) we calculate the 
spin and orbital components of the ordered moment to 
be 2.69 hb and 1.12 fiB, respectively. The unquenched 
orbital component has also been shown to be impor- 
tant for achieving a quantitative understanding of the 
bulk susceptibility— The total ordered moment in the 
model of 3.8/xb is rather larger than the experimentally- 
determined ordered moment of 2.9 /xb (Ref. @). It is not 
clear why these values should differ, but if confirmed then 
this difference could be the reason why the calculated 
spectrum has a higher overall intensity than the mea- 
sured spectrum. 

Perhaps the most surprising result to emerge from 
the analysis of the spin excitation spectrum is the vir- 
tual absence of the 90deg Co 2+ -Co 3+ -Co 2+ exchange 
couplings Ji and J2 relative to the 180 deg coupling J. 
This is dramatically counter to the simplest estimate 
Ji ~ J 2 ~ 2 J based on exchange paths involving a 
single Co 3+ bonding orbital^ and implies that frustra- 
tion effects on the magnetic order are very small. One 
consequence is that the magnetic structure can be re- 
garded as two interpenetrating, square-lattice, antifer- 
romagnets with spacing 2a. In the complete absence 
of Ji and J2 the moment directions for the two inter- 
penetrating antiferromagnets would be unrelated to one 
another, but quantum fluctuations are expected to sta- 
bilize collinear order of the global magnetic system via 
the order-by-disorder mechanism.— A study of the spin 
wave spectrum of stripe-ordered La2-a;Sr :r Ni04 similarly 
concluded that the 180 deg Ni 2+ -Ni 2+ exchange was sig- 
nificantly larger than the 90 deg coupling. 42 Given the 
importance of exchange for the stability of spin- and 
charge-ordered ground states it would be of interest to 
seek an understanding the exchange interactions in these 
systems in terms of the bonding orbitals involved. 

There remain some details of the magnetic order to 
be understood. In Fig. [5^a) we reported how the ratio 
Jiio/Aio increases from approximately 0.4 to almost 1.0 
in the temperature range 30— 50 K, and in Fig. [4] we see 
that the start of this increase coincides with a kink in 
the magnetization below which the field-cooled and zero- 
field-cooled magnetization separate. To these observa- 
tions we add a recent finding from a muon-spin rotation 
(/iSR) study^ 3 - that the kink coincides with the temper- 
ature at which magnetic ordering sets in on the muon 
timescale (~ 10~ 6 s). Since /iSR and bulk magnetization 
probe much longer timescales than neutron diffraction 
(~ 10~ 12 s) we can take 31 K to be the static magnetic or- 
dering temperature. Between 31 K and ~ 60 K magnetic 
Bragg peaks are still observed by neutron diffraction and 
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so in this temperature range the magnetic order is not 
static but fluctuates on a timescale between ~ 10 -6 and 
~ l(T 12 s. 

For the assumed collinear magnetic structure there are 
several possible ways to explain the behavior of iiio/JiTo- 
One is in terms of a canting of the moments. If the mo- 
ments point along the a or b direction then iiio/-Tiio = 1, 
but if they make an angle greater than 45deg to [110] 
then iiio/Alo < 1 [note that we are referring here to the 
domain in which the modulation of the magnetic struc- 
ture is along [110] — see Fig. HJb)]. The observed low 
temperature value of Jiio/-fiio = 0.4 corresponds to a 
turn-angle of 12deg. According to this model, therefore, 
the moments point along the a (or b) direction when they 
first start to order, but gradually turn towards the [110] 
direction (i.e. away from the modulation direction) on 
cooling. By ~ 30 K the moments have turned through 
an angle of about 12deg, i.e. in Fig. |ljb) <f> = — 12dcg 
or +102 deg, and they remain at this angle at lower tem- 
peratures. Although consistent with the data, this model 
requires an explanation for what causes the rotation of 
the ordered moments and for why they fix on a canting 
angle of 12 deg from the Co-0 bond directions at low 
temperatures. 

Another possible interpretation of the Jno/Jiio data 
is in terms of a change in population of different inequiv- 
alent spin domains. For example, if the easy magnetic 
direction within the plane were at 45 deg to the Co-0 
bonds, i.e. 4> — ±45 deg in Fig. [IJb), then in the or- 
dered phase the moments could in principle point either 
parallel or perpendicular to the in-plane modulation di- 
rection of the structure. These longitudinal and trans- 
verse structures are not related by symmetry so their 
energies will in general differ. If the energy of the trans- 
verse structure were the lower of the two then the trans- 
verse structure would be present in greater proportion 
at the temperature of ~ 30 K when static order sets in, 
resulting in -/no/Aio < !• With increasing temperature 
above ~ 30 K thermal fluctuations would tend to equalize 
the populations of the two domains, consistent with the 
observed increase in Iuq/I^q. The difference between 
the field-cooled and zero-field-cooled magnetization be- 
low 30 K may be caused by the effect of the magnetic 
field favoring one domain over the other. 

Keeping these models for the magnetic order in mind 
let us now consider the data on the orthogonal compo- 
nents of the in-plane magnetic fluctuations at 2 meV and 
4meV ( Section [III Cp . Assuming transverse fluctuations 
of the moments about an angle <fi = —12 deg with re- 
spect to the Co-0 bonds we expect to measure a ratio 
S 110 /S 1W = cot 2 (45° + 4>) = 2.4, whereas experimentally 
we find ratios of 1.19 ±0.07 at 2meV and 1.00 ±0.05 at 
4meV corresponding to an angle of at most 1 or 2 deg 
away from the Co-0 bond direction. Because the scatter- 
ing is proportional to the square of the orthogonal com- 
ponents of the moments the model with unequal domain 
populations predicts the same ratio S 110 / S 110 = 2.4. 
Both models, therefore, are apparently at odds with the 



data. 

However, the elastic and inelastic scattering data can 
be reconciled through consideration of the particular 
form of the magnetic order which, as discussed above, can 
be regarded as two interpenetrating, square-lattice anti- 
ferromagnets with spacing 2a. In isolation, each of these 
antifcrromagnets would have a magnetic dispersion with 
four-fold symmetry so that, for example, there would 
be equivalent minima in the dispersion at each of the 
in-plane wavevectors ±(0.25,0.25), and ±(0.25,-0.25). 
Weak coupling between the two antiferromagnets locks 
them into a single structure with two domains each hav- 
ing a two- fold pattern of magnetic Bragg peaks, but the 
four-fold symmetry will remain in the excitation spec- 
trum above a crossover energy related to the coupling 
between the two antiferromagnets. This has two conse- 
quences for the data here. One is that after averaging 
the excitation spectra from the two wavevector domains 
the in-plane fluctuations are isotropic above the crossover 
energy. Therefore the observation of S 110 / S lw ~ 1 
at 2meV and 4meV suggests that the crossover energy 
is below 2 meV and is further evidence that the mag- 
netic order comprises two weakly-coupled, interpenetrat- 
ing antiferromagnets. The other consequence is that the 
magnon dispersion surface above the crossover energy is 
the same for both wavevector domains and corresponds 
to a square-lattice collinear antiferromagnet with spac- 
ing 2a. This means that domain-averaging and the fact 
that we do not have a unique model for the magnetic or- 
der has no bearing on the interpretation of the excitation 
spectra. 

Finally, we note that we have not found any evi- 
dence for magnetic degrees of freedom associated with 
the Co 3+ site in the magnetic spectrum probed here up 
to ~ 50meV. However, we cannot rule out the possi- 
bility of a small Van-Vleck moment on the Co 3+ sites 
induced by coupling to the Co 2+ moments. If this were 
the case then such a coupling might be able to cause a 
spin-canting transition, and this might be another possi- 
ble explanation for the -/no/Aio data. In the absence of 
experimental evidence to test this possibility we do not 
attempt to speculate further, but since other spin-charge 
ordered systems exhibit similar changes in /no/Aio with 
temperature to that found her o 10 ' 29 ! 30 ' 31 it would be of 
interest to examine this behavior more closely. 



VI. CONCLUSIONS 

We have gained a rather complete understanding of 
the nature of the magnetic excitations in what is a text- 
book charge-ordered, two-dimensional antiferromagnet. 
The magnetic order is stabilized essentially by a single 
exchange interaction acting along a straight-line path be- 
tween the charge-ordered Co 2+ sites. We find no evidence 
for active magnetic degrees of freedom on the Co 3+ sites. 
Open questions include what is precise nature of the mag- 
netic order and how to explain the exchange in terms of 
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the bonding. 



B. Many- level spin— orbital model 
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APPENDIX: GENERALIZED LINEAR 
SPIN- WAVE THEORY 

A. Effective spin— | model for ground state doublet 

Quantization of Eq. ^ by means of the Holstein- 
Primakoff transformation, followed by diagonalization of 
the resulting Hamiltonian by the standard method^ 
leads to two non-degenerate branches with in-plane dis- 
persion 



S^i,2(Q) = y (^4q ± Bq) 2 - £>q, (4) 
where u)\ (0J2) corresponds to the upper (lower) sign, and 

Aq = 25{2J(l + e) + J1-J2 + J 2 cos[Q- (a- b)]} 

B Q = S{J6cos(2Q-a) + J<5 cos(2Q • b)}, 

Dq = 2S{ J(l - 5/2)[cos(2Q • a) + cos(2Q • b)] 

+Jicos[Q ■ (a + b)]}. (5) 

With the spins aligned along x, only the transverse corre- 
lations yy and zz contribute to the linear spin-wave cross 
section, Eq. ([1]). The corresponding response functions 
(per Lai.5Sro.5Co04 formula unit) for magnon creation 
are given by 



S zz {Q,u) 



qIsa 



We diagonalizc (J3]) in two steps. First, we diagonal- 
ize the single-ion terms, i.e. the crystal field and spin- 
orbit terms, plus the molecular field part of the exchange 
energy. This produces a set of self-consistent single-ion 
energy levels and wavefunctions for each site. The sec- 
ond step is to write the residual part of the exchange 
interaction in terms of pseudo-boson raising and lower- 
ing operators for the single-ion states, retaining terms up 
to quadratic order. The resulting Hamiltonian is bilinear 
in the pseudo-boson operators and can be diagonalized 
by the standard procedure. 

To be specific, let us assume the antiferromagnetic or- 
der in Lai 5Sro.5Co04 to be composed of two sublattices 
A and B, with the ^4-sublattice moments along +x and 
the i3-sublattice moments along — x. The Hamiltonian is 
the sum of a single-ion part Hi and a two- ion part H.2- 
The single-ion Hamiltonian for one unit cell is given by 



rti = rtj + rt 1 



where 



njA njA 

n 1 - n C f 



H s . 

Hi 



H 



ml' 



(7) 



(8) 



In the first of these equations H^, T~t^ and H^ f are the 
crystal field, spin-orbit and molecular field interactions 
for the A site. The latter is given by 



H mf 



{A B } 



E 



(9) 



Q 



Q 



D 



Q 



4 fiwi(Q) 
g 2 z SA Q -B Q -D Q 



{n(u) + l}6[w-wi(Q)], 

{n(«) + l}*[w-(^(Q)], 
(G) 



A and B are interchanged for the B site terms in ©. 
Aa and represent the displacements from an A site 
to other A and B sites, and Ja a and Ja b are the corre- 
sponding exchange parameters. In practice, we restrict 
the model to nearest-neighboring sites, so the summa- 
tions over the B- and ^4-site neighbors in Eq. ([9]) amount 
to 4J + 2Ji and 2J 2 , respectively (see Fig.[lJ. Note that 

(s B )--(s^>. 

The mean-field equations ([5]) are solved self- 
consistently by an iterative process until the values of 
(S A ) and — (S B ) converge to an acceptable level of pre- 
cision (in our case one part in 10 6 ). This gives the set of 
single-ion energy levels = e B = e„, where n takes val- 
ues from (ground state) to (2L + 1) x (25+1) - 1 = 27. 
From the corresponding single- ion wavefunctions |n), the 
matrix elements for spin 



where g y and g z are in-plane and out-of-plane ^-factors 
for the effective spin-| ground-state doublet of the Co 2+ 
ion and n(u) is the boson occupation number. We neglect 
the corrections which are sometimes applied to the spin- 
wave dispersion and response functions to account for 
zero-point fluctuations. 



S n 'n = (n'\S\n) 
and for the total magnetic moment 

M n , n = -<n'|L + 2S|n) 



(10) 



(11) 



can be calculated for the A and the B sites. 

We now consider the two-ion part of the Hamiltonian 
which describes the residual exchange interactions. By 
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symmetry, the interaction energy is the same for the two 
sites, so we can write it as twice the energy for the A site, 

H 2 = ]T J A S A ■ S A - 2S A ■ U A { 

{A} 

= E j a (s a - (s A )) ■ (s A - <s A )) - (s A ) ■ n A , . 

{A} 

(12) 

The summation is over all sites connected to the A site by 
non-zero exchange interactions. The molecular field part 
of the exchange interaction is subtracted because this is 
included in the single-ion terms — see Eqs. (jSJ. 

To quantize the Hamiltonian we introduce pseudo- 
boson raising and lowering operators which convert the 
ground state into the excited states, and vice- versa. For 
the A site, 

4|0) = \n) and a n \n) = |0) . (13) 

Operators and b n are defined similarly for the B site. 
If the temperature is sufficiently low that the equilibrium 
population of the ground state is close to one then to 
a good approximation these operators satisfy the Bose 
commutation relations^ 



Operators on different sites commute. The single-ion 
Hamiltonian can now be written as 

H\=} j e n (ala n + b\b n ) . (15) 

n>0 

For the two-ion Hamiltonian we start with the follow- 
ing identity for the spin operator on the A site, 

S = Sqq + y^[SnQQ^ + So n a n + (S nn — Soo)al l a n ] 

n>0 

+ E E S n / n a^,a n . (16) 

At low temperatures Soo ~ (S), and we need only retain 
the linear terms in the operators if we are to neglect 
higher-than-quadratic terms in the Hamiltonian. With 
these approximations Eq. (|16[) simplifies to 

S - (S) = £ Snoai + S 0n a n , (17) 

n>0 



[a n ,a n ,] = [b n ,bl,] = 8 nn i . (14) and the two-ion Hamiltonian (|T2|) becomes 



H 2 



(S A ) ■ H 4 



f t E E J Aa[( s «o • So„<)4 a n',A A + (S^ • S A )al l a f n , AA +h.c] 

n,n'>0 A.a 



n,n'>0 A B 

I 

The Fourier transform operators are defined by be written in the form 



(18) 



wa = -^E ex p[*Q' ( m+ A )i 

* q 

4+a = ^E exp [~ iQ ' ( m+ A 



)] < , (19) 



where TV is the total number of A sites (or B sites). In 
Eq. (UHJ) we explicitly show the position vectors for the 
operators: m is position vector for the A site in the m th 
unit cell, so for example al n+Ai is the raising operator for 
the A site that is displaced from m by A . The definitions 
of the Fourier transform operators for the B site are the 
same as those for the A site except b replaces a. After 
substitution of the expressions in Eqs. ([TO]) into (fT5")) and 
(I18p and summation over m the total Hamiltonian can 



TL = Ho 



^E E 

Q n,n'>0 



(20) 



where 7io contains the constant terms, ~X. n q is the row 

matrix (a' n q, b' n q, a„._Q, 6„.-q) for the excited level n, 
X^q is the column matrix containing the Hermitian ad- 
joint operators, and 



'nn' ,Q 



^ A nn ' 

B* 

c: 



Q 

n'n.Q 



',Q C nn '.Q D nn i t Q 

',Q fn'n Cnn',Q 

',-Q Ai'n,Q B* m , _q 

Cnn',Q Bnn',-Q A n >n,Q / 



(21) 
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The coefficients in the matrix are 

A nn >,Q — e n 5 nn i + (S^ • Sg n ,)7q , 
B nn ',Q = (S^ • S^„,)7q , 

C 7l n',Q = (S n0 ' S„/q)7q , 

Azn',Q = (S^ • S^, )7q , (22) 

where 

7q = J]] Ja a exp(iQ • A A ) , 

7q = J ^b ex P (iQ • A B ) . (23) 

In general, 7^q = 7q because all A sites are equivalent 
on the magnetic lattice. For the present system, 

7^ = 2 J 2 cos[Q ■ (a - b)] , 

7^ = 2Jcos(2Q • a) + 2Jcos(2Q • b) 

+2 J a cos[Q • (a + b)] . (24) 

The Hamiltonian (|20|) can now be diagonalized by the 
standard method.— There are a total of (2L + 1) x 
(25 + 1) — 1 = 27 single-ion excited levels for Co 2+ in 
Lai.5Sr .5CoO4 giving a total of 54 distinct modes in the 
magnetic spectrum, two modes for each single-ion excited 
level n. In our case we diagonalized the full Hamiltonian 
which is represented by a 108 x 108 matrix (since each 
mode appears twice in the Hamiltonian). 

To evaluate the neutron scattering cross section we em- 
ploy the general form for the response function that takes 



into account orbital as well as spin magnetization. For 
the creation of one magnon in the \rij) mode (j = 1,2) 
from the fully ordered ground state via aa correlations 
the response function (per Lai.sSro.sCoC^ formula unit) 
is given by 

S aa (Q,u;) = 1 | (nj\M a (Q) |0)| 2 8[w - u nj (Q)] . 

(25) 

By replacing S by M in Eq. (117|) and summing over the 
A and B sites we obtain the following expression for the 
operator representing the Fourier transform of the mag- 
netization, 

M°(Q) = J2 [ M,)Xq + WJ V-Q 

n>0 

+W) B ( Q + (C)\-Q]- 

(26) 

The a n q, a ni _Q, etc., operators are expressed as a linear 
combination of creation and destruction operators for the 
magnon modes via the Bogoliubov transformation ma- 
trix. From the coefficients of the creation operator for 
the \rij) mode one can calculate the matrix element in 
(|2"6")) . and hence the response function Eq. (|2"B"|) and scat- 
tering cross section Eq. ([T]) for magnon creation in this 
mode. The resulting scattering intensity and response 
functions for the lowest six modes of Lai.sSro.sCoC^ are 
shown in Figs. [7] and[H] 
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